function F = dfun(x)
F = [2*x(1) - 2 ,-1;2*x(1) , 8*x(2)];
end
